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Introduction: To explain the equatorial valley 
networks on Mars, Jakosky and Carr [1] suggested that 
water ice now stored in the north polar region would 
be mobilized at high obliquity and precipitate out at 
low latitudes. Extrapolating the present day latitudinal 
distribution of water vapor to high obliquity condi- 
tions, and noting that the low latitude atmosphere 
would be saturated, they predicted substantial surface 
ice deposits would accumulate in the tropics at such 
times. 

The first general circulation model simulations to 
verify this prediction were reported by Haberle et al. 
[2] who found that while ice can accumulate at low 
latitudes at high obliquity, it is distributed regionally 
depending on orbital conditions. Forget [3], 
Richardson and Wilson [4], and Mischna et al. [5], 
subsequently obtained similar results with independent 
models. Thus, obliquity driven climate change may 
help explain the many tropical landforms thought to be 
sculpted by water in one form or another (see, for ex- 
ample, refs [6], [7], and [8]). 

While low latitude ice accumulations at high obliq- 
uity appears to be a robust result, the major challenge 
now facing models is predicting ice accumulations in 
the same places where the geological evidence sug- 
gests it occurred. This will depend not only on orbital 
conditions, but also on what physical processes the 
models include in the hydrological cycle. For example, 
none of the models mentioned above include the ra- 
diative effects of water vapor or clouds, yet both are 
expected to be in abundance at high obliquity. And 
none of the models has a very realistic cloud micro- 
physics scheme, which can have a significant effect on 
how clouds affect the planet's radiation balance. 

Here we extend these early modeling results by in- 
cluding a more sophisticated cloud microphysics pack- 
age, as well as the radiative effects of water vapor and 
clouds. 

Model description: We use the NASA/Ames C- 
grid Mars general circulation model with an updated 
radiation code and cloud microphysics scheme. To 
speed up the simulations, we run the model at fairly 
coarse resolution (7.5° latitude x 22.5° longitude). Fu- 
ture efforts will examine the effect of resolution on the 
results. 


Radiation Code Fluxes and heating rates are cal- 
culated from a radiation code based on the two-stream 
solution to radiative transfer that fully accounts for 
multiple scattering in the presence of gaseous absorp- 
tion. The model has 12 spectral intervals. Dust and 
water ice scattering properties are included. For dust, 
we use the Ockert-Bell [9] values in the visible, and 
Forget [10] values in the infrared. For ice, we can ei- 
ther compute them online as the cloud evolves, or we 
can specify them. Gaseous opacities for water vapor 
and CO 2 are calculated from correlated k-distributions 
taken from full line-by-line models. 

Cloud Microphysics Our cloud scheme is based on 
a moment/order scheme in which the mass mixing ra- 
tio and number density of the cloud ensemble are the 
advected species. From these we obtain a mean particle 
size and an estimate or the particle size distribution 
(assuming a variance) which we then divide into 8 
bins. Cloud microphysics is performed in each of these 
bins and includes nucleation, condensation, and gravi- 
tational settling. Dust is treated as a tracer and serves 
as condensation nuclei. The altered size distribution is 
then converted back into a mean size, a mixing ratio, 
and a particle number density. 

Results: We have conducted simulations for a va- 
riety of different obliquities, all at present solar lumi- 
nosity. In each case the model is spun up from dry ini- 
tial conditions with a residual ice cap at the north pole. 
After several years, depending on obliquity, the atmos- 
phere equilibrates and repeats from year-to-year. A 
sample result for the 6(T obliquity simulation, without 
the radiative effects of clouds or water vapor, is shown 
in Fig. 1. The top panel in Fig. 1 is the zonally aver- 
aged column water vapor as a function of time for 7 
Mars years. The middle and bottom panels are similar, 
but for cloud mass and surface ice, respectively. 

Water ice subliming from the north residual cap 
during summer is rapidly transported southward. 
Clouds form in low northern latitudes and ice precipi- 
tates to the surface. The remainder is transported into 
the southern hemisphere and condenses onto the south 
seasonal CO 2 ice cap which extends almost to the 
equator at the solstice. When the south cap retreats, 
water is released into the atmosphere where some pre- 
cipitates back to the surface and the remainder is trans- 
ported north. Again clouds form in the low latitudes 



Third Mars Polar Science Conference (2003) 


8060.pdf 


and ice precipitates to the surface. At equilibrium, 
thousands of precipitable microns of water vapor ap- 
pear in the summer polar regions. There is more water 
in the south than the north because the south cap is a 
better trap for water, and because the Southern Hemi- 
sphere is warmer during summer than in the north. 
Cloud abundances also reach the thousand precipitable 
micron mark with model predicted particle sizes in the 
20-30 micron range. These particles are much bigger, 
and subsequently fall out faster, than those for present 
obliquity. 

Eventually, permanent deposits form (i.e., ice re- 
mains on the ground all year long) in the low latitudes 
of each hemisphere. These deposits are concentrated 
along the northern flanks of the Tharsis region and to 
the northeast of the Hellas basin. Topography plays a 
key role on where the deposits form through its influ- 
ence on the circulation. The deposits do not necessarily 
form in locations where the mean annual surface tem- 
peratures are a minimum. They form where the satura- 
tion state of the atmosphere is highest. This, in turn, is 
influenced not only by the thermal structure of the at- 
mosphere, but also by the transport characteristics of 
the atmosphere. 

Simulations which include the radiative effects of 
water vapor show similar results, but with (a) an in- 
crease in the amount of surface ice, (b) a slight shift in 
the location of the deposits, (c) a cooler and cloudier 
atmosphere, and (d) slightly warmer surface tempera- 
tures. We are presently undertaking simulations with 
the radiative effects of clouds included and will report 
the results at the meeting. However, off line 1-D 
simulations using the predicted cloud abundances indi- 
cate they will have a much greater influence on the 
results than water vapor alone. Their abundances (~ 
1000 pr-pm), particle sizes (20-30 pm), widespread 
occurrence, and impact on the solar and infrared radia- 
tion fluxes give clouds a much greater role in deter- 
mining the climate at high obliquity than for present 
day conditions. 

Conclusions: Mars has a natural mechanism for 
experiencing significant climate change and redistrib- 
uting surface ice. Obliquity changes alone are quite 
capable of moving ice into low latitudes and may pro- 
vide an explanation for the many geological landforms 
that strongly indicate recent climate change. 
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